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Abstract. Many embedded star clusters do not evolve into long-lived bound clusters. The most popular explanation for this 
"infant mortality" of young (few Myrs) clusters is the expulsion of natal gas by stellar winds and supernovae, which perturbs 
the clusters' potential and leaves up to 90% of them unbound. A cluster disruption model has recently been proposed in which 
this mass-independent disruption of clusters proceeds for another Gyr after gas expulsion. In this scenario, the survival chances 
of massive clusters are much smaller than in the traditional mass-dependent disruption models. The most common way to 
study cluster disruption is to use the cluster age distribution, which, however, can be heavily affected by incompleteness. To 
avoid this pitfall we introduce a new method of studying cluster disruption based on size-of-sample effects, namely the relation 
between the most massive cluster, M max , and the age range sampled. Assuming that clusters are stochastically sampled from a 
power-law cluster initial mass function, with index -2 and that the cluster formation rate is constant, M max scales with the age 
range sampled, such that the slope in a log(M max ) vs. log(age) plot is equal to unity. This slope decreases if mass-independent 
disruption is included. For 90% mass-independent cluster disruption per age dex, the predicted slope is zero. For the solar 
neighbourhood, SMC, LMC, M33, and M83, based on ages and masses taken from the literature, we find slopes consistent with 
the expected size-of-sample correlations for the first 100 Myr, hence ruling out the 90% mass-independent cluster disruption 
scenario. For M51, however, the increase of log(M max ) with log(age) is slightly shallower and for the Antennae galaxies it is 
flat. This simple method shows that the formation and/or disruption of clusters in the Antennae must have been very different 
from that of the other galaxies studied here, so it should not be taken as a representative case. 

Key words, galaxies: star cluster - galaxies: evolution 
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1. Introduction 

In this study we propose an alternative method of studying the 
nature of cluster disruption, by using only a small number of 
massive clusters and a few simple assumptions. This study is 
motivated by the recent debate on the duration of the "infant 
mortality" phase of clusters. 

The term infant mortality was coined by Lada & 



Lada d2003l LL03), who noticed from a comparison of the 
number of embedded clusters per unit time (dN/dt) to the 
dN/dt of optically detected clusters in the solar neighbour- 
hood, that up to 90% of the embedded clusters do not sur- 
vive the gas expulsion phase. This rapid expulsion of gas is 
driven by stellar winds, ionisation and supernovae of early type 
stars. Due to the expanding gas shells the binding energy of 
the initial system of stars and gas is reduced and the stars that 
have remained in place (assuming instantaneous gas expulsion) 
suddenly have velocities higher than the loca l escape velocity 
dTutukovll 19781: lHiUj|l980ULada et alJ ll984t lGoodwinlll997: 
Gever & Burkertl 200 if Kroupal 2005 ). Under the assumption 



of a constant, or mass-independent, star formation efficiency, 
the fraction of clusters that becomes unbound (infant mortal- 
ity rate), or the fraction of mass lost from each cluster ("infant 
weight loss") is independent of the mass of the embedded clus- 
ter. Theory predicts that the effects of gas expulsion s hould be 
largely over within 20 Myr ( Goodwin & Bastian1 l2006). 

More recently, the infant mortality scenario has also been 
used to explain the steep drop in dN/dt around ~ 10-20 Myr 
of clusters in M51. Around 70% of the clusters in M51 do 
not survive past 20 Myr, roughly inde pendent of cluster mass 
l|Bastian etal J 120051: iGieles et al]l2005l) . It is noteworthy that 
the study of M51 clusters was based on optically detected clus- 
ters only, while the study of the solar neighbourhood (LL03) 
used a comparison between young embedded and older opti- 
cally visible clusters to determine the infant mortality fraction. 
The dN/dt distribution of optically detected cluster s in the solar 
neigh bourhood is nearly flat for ages < 1 00 My fl ( La mers et al 



2005a), i.e. very different than that of the M51 clusters. 



Send offprint requests to: mgieles@eso.org 



1 For a sample limited to clusters within 600 pc from the sun and 
with masses > 100 M a . When plotting dN/dt of all clusters, as in 
Whitmore et al. ( 2007), it declines due to incompleteness effects. 
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Fall et al.1 d2005l) find that the dN/dt distribution of a 
mass limited cluster sample of the Antennae galaxies declines 
roughly as t~ l up to ~ 1 Gyr. They explain this decline by 
infant mortality, removing 90% each age dex independent of 
cluster mass over the full age range of the dN/dt distribution. 
Since this time-scale is two orders of magnitudes longer than 
the time-scale involved in the original definition of infant mor- 
tality by LL03, we will refer to this disruption model as "mass- 
independent dissolution" (MID). 

Based on the Magel lanic Cloud Photometric Survey, 



Rafelski & Zaritskvl d2005l) presented a set of structural param- 
eters of clusters in the Small Magellanic Cloud (SMC). Their 
version of th e dN / dt distribution w as declining and this was in- 
terpreted by IChandar et al. (2006) as MID removing the same 
fraction (90%) of clusters each age dex as in the Antennae dur- 
ing the first 1 Gyr. 



(2006) 



is surprising 



The interpretation of IChandar et al. 
for two reasons: a) it is in disagreement with earlier stud- 
ies on the dN/dt distributio n of SMC clusters, who found a 



Chiosi et alj 120061) 



flat dN/dt up to ~ 1 Gyr dHodgdll987 : 
and b) they are the first to suggest a disruption scenario in 
which the life-time of star clusters during the first Gyr does 
not depend on their mass or local environment. This con- 
tradicts the existing theoretical understanding of cluster dis- 
ruption. Clusters in the SMC should survive much longer 
than clusters with similar masses in the Antennae galax- 
ies, due to the much lower tidal fiel d strength and molec 



ular c loud density in the SMC (e.g. ISpitzerl 119581 IWielenl 

"It 



1 98 81 1 Vesperini & Heggid 1 1 9971: iBaumgardt & Makinoll2003l 
Gieles et a l. 2006c). From comparisons between the age and 
mass distributions of clusters in different gala xies this sce- 
nario is also sup 



nano is also supported by observations (e.g. lElson & Fall 
19851: iHodgdl 19871: iBoutloukos & Lamersll2003l [Tamers et al 
20053). The decline of the dN/dt distribution of SMC clusters 



reported by Chandar et al. ( 20061) is most probably the result 
of detection incompleteness since they derive the slope of the 
dN/dt distribution from a fit to t he ful l cluster sample, without 
making a mass cut. Gieles et al. ( 20071) showed that a mass lim- 
ited sub-sample has a flat dN/dt up to 1 Gyr. This is because at 
older ages clusters are intrinsically fainter, causing the number 
of observed clusters in the full (luminosity limited) sample to 
drop with increasing age. 

The interpretation of results based on dN/dt distributions 
will always be heavily dependent on how incompleteness is de- 
termined or corrected for. To remedy these shortcomings we 
have developed a method that can serve as an in dependent 



check of the universa l MID scen ario proposed b y Fall et al 



d2005l) . IChandar et all d2006l) and IWhitmore et alj d2007l) and 
the traditional disruption scenario in which the disruption time 
depends on mass. 

Our new method is based on a few very elementary as- 
sumptions and needs only a handful of massive clusters at var- 
ious ages, i.e. well above the detection limit, avoiding prob- 
lems with incompleteness of faint clusters at old ages. An ad- 
ditional advantage of our method is that age dating of massive 
clusters is more accurate than for clusters with lower masses 
since the stellar IMF is well populated, so stochastic fluctua- 
tions in the cluster colours due to IMF sampling are small. The 



meth od is similar, but unde r different assumptions, to that used 
by lMaschberger & Kroupa (2007) who used the most massive 
cluster as a function age to derive the star formati on history of 



galaxi es. It is fundamentally similar to the study of Hun ter et al 



(2003) who used the most massive cluster per logarithmic age 
bin to constrain the initial mass function of clusters. 

Using the size-of-sample effect of cluster populations, in 
particular the most massive cluster per logarithmic time inter- 
val, we will address the following points in this paper: 

- Is cluster disruption independent of cluster mass and envi- 
ronment? 

- How does the observed relation between the most cluster 
found and the linear time interval probed depend on cluster 
disruption and the cluster formation history of the galaxy? 

- What is the relative fraction of stellar mass in bound and 
long lived (> lOMyr) clusters? 

The proposed method is meant to compliment, and provide 
an independent check, of existing methods, which are normally 
based on the observed age distribution of full cluster popula- 
tions. 

In §|2] we present the statistical considerations of the size- 
of-sample effect and compare the expected behavior to ob- 
served cluster populations in § [3] A discussion and our con- 
clusions are presented in §|4]and §|5] respectively. 

2. The maximum cluster mass from statistical 
considerations: size-of-sample effect 

2.1. The mass of the most massive cluster 

We assume that masses of clusters follow from random sam- 
pling of a power-law cluster initial mass function (CIMF): 



dN 

MM) = =AM-'\ 

VV ' dM 



(1) 



with or 2 dZhang&Falj|l999tlBiketal.ipi 
20031: Ide Griis et alJl2003HGieles et al.ll2 006a). 



2003; Whitmore 



For a < 2 the total mass diverges to infinity when integrat- 
ing to infinity, so a maximum mass has to be chosen, which can 
be interpreted as a physical maximum above which no clusters 
can exist. We will refer to this mass as M up (i.e. the upper limit) 
and to the most massive cluster observed, i.e. the most massive 
cluster actually formed, as M max . 

The value of M max depends on the constant A in Eq. Q]and 
can be found by integrating <p(M) from M max to M up and setting 
this equal to 1. For a > 1 and M up >> M max , which is probably 
true for most galaxies, this results in A = (a - lJM^ 1 . For 
a = 2 this reduces to A — M max . The relation between M up 
and M max depends on the number of clusters formed (AO which 
is proportional to A. The larger N, the closer the statistically 
probable M max will be to M up . We can relate M max to N as 



N 



JM mh 

max 



A M~ a dM 



,a > 1 



(2) 
(3) 
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Mm 

M„ 



(4) 



where in the last steps we again assumed M up /M max >> 1. For 
a constant M m ; n we see from Eq.[3]that M max scales with N as 



rl/(<*-l) 



(5) 



The sca ling of M max with N for clusters (Eq. [5]) is analogous 
to what lOev & Clarke (2005) found for stars. Assuming a con- 
stant cluster formation rate (CFR), dA7df=constant, the num- 
ber of clusters in each equal logarithmic time interval, Nbm, 
increases linearly with age, since 



N h 



dN 
dln(f) 



dN 
~d7 



(6) 



and so A r bi n (0 K ?■ The same holds for dAVdlog(f) apart from 
an additional constant ln(10). Therefore M max in logarithmic 
age bins scales as 



Mrr 



oc N l J {a - l) oc t lKa - X \ 
bin 



(7) 



The above derivation was based on a constant CFR and no dis- 



ruption. iHunter et al.l d2003l) already noticed that a power-law 
relation for the CFR with age (CFRoc t n , with 77 negative for a 
CFR that was lower in the past) would change the relation for 



M m:ix to 



(8) 



The scaling of M max with age appears since Nb m increases 
in equal sized logarithim age bins. O ur assumptions are fun- 
damentally different from those of Maschberger & Kroupa 
(2007), since in their model M max is set by the star formation 
rate and a time-scale of formation of an entire cluster popula- 
tion of lOMyr. Though they do randomly sample M max from 
a probability density function, allowing some influence of the 
size of the sample, they on average predict a constant M max 
in logarithmic age bins for a constant CFR, rather than the in- 
crease predicted by Eq. [7] Also, they assume a cluster disrup- 
tion law and solve for the cluster formation history. We, on the 
other hand, take the cluster formation rate as constant and test 
a theory of cluster disruption. 

For the moment we consider the simplified scenario of a 
constant CFR and no mass-dependent disruption (Eq.lTJ. In the 
next section (§ \2.2\ we will add the MID scenario to this. In 
§ |4] we discuss the effects of mass-dependent disruption and 
variations in the cluster formation history. 

2.2. The effect of mass-independent dissolution on 

M max 



We consider t he sc enari o proposed by iFallet alJ (12005b . 
Chandar et ail J2006I) and IWhitmore et al.l d2007b . in which 
each age dex a fixed fraction of the number of clusters gets de- 
stroyed, independent of the cluster mass. We refer to this frac- 
tion as /mid - They argue that 90% of the clusters are destroyed 
each age dex, so /mid = 0.9. This reduction in number results 
in an expression for the remaining clusters as a function of time 



0) 
< 



Q. 

O 




0.5 



0.0 

0.0 0.2 0.4 0.6 0.8 1.0 

Fraction destroyed each age dex (f M!D ) 

Fig. 1. Predicted slopes for log(M max ) vs. log(age) for different 
MID fractions (/mid) and three different indices of the CIMF 
(a) (Eq.HU). 



of the form AN I At oc t A , with A = log(l - /mid), which is r 1 
for /mid = 0.9. For/ M iD=0, 0.5 and 0.8 we find A _ 0, -0.3 and 
-0/7, in agreement with the values quoted by Whi tmore et al. 
(2007). From this relation for dN/dt and Eq.[6]we find that the 
number of clusters per logarithmic age bin, Nu n , depends on 
/mid as 



AW <* 



dN 
dlog(f) 



oc t 



l+log(l-/ M i D ) 



(9) 



We now substitute this expression in Eq. [5] for N to find the 
trend of M max with age when MID is included: 



oc t 



(l+log(l-/„ ID ))/(a-l) 



(10) 



which gives the result as in Eq. [8] for no disruption (/mid = 
0) and a constant CFR (77 = 0). From a comparison between 
Eq.[8]and Eq.[10]it is directly visible that there is a degeneracy 
between a CFR that has been increasing (rj < 0) and MID. 
MID reduces the number of clusters with increasing age. Due 
to the mass-independent nature of this disruption model, it is 
impossible to distinguish between an increasing CFR and MID. 
In a plot of log(M max ) vs. log(age) we thus expect a slope 



slope 



1 +log(l - /mid) 



a- 1 



+ 1 





for /mid 
for /mid 



: 0, a = 2, 
0.9 and all a. 



(11) 
(12) 
(13) 



In Fig. Q] we show the predicted slopes following from 
Eq.QTjfor three different values for a. For the universa l /mid 



0.9 scenario, which was proposed by Whitmore et al. 1 200% 



the predicted slope is for all values of a. For lower values 
of /mid, we expect noticeably different slopes than +1, provid- 
ing an independent strong constraint of the acceptable values 
of /mid- 
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3. Comparison to the observations 

3.1. Description of the data used 

We collect cluster ages and masses in seven different galax- 
ies from the literature: the Milky Way (solar neighbourhood), 
SMC, LMC, M33, M83, M51 and the Antennae galaxies. Here 
we briefly mention the origin of the data and we refer the reader 
to these papers for details on the data reduction and age fitting 
techniques. 

For the clusters in t he solar neighbour h ood w e use ages 
from the catalogue of Kharchenko et al.l d2005l) and cor- 
responding masses derived and kindly provided to us by 
Lamers et al.l d20Q5al) . The mass estimates are derived from the 
number of stars with high membership probability and an ex- 
trapolation of the stellar initial mass function. We limit our- 
selves to the 209 clusters in the catalogue that are within a dis- 
tance of 1 kpc from the sun, for which the mass estimates are 
believed to be accurate and the sample is not severely affected 
by distance incompleteness. 

For the SMC and the LMC we use the results of 
who kindly provided us with a table 
of 191 SMC clusters and 748 



Hunter et al 



(120031) 

of ages and luminosities 
LMC clusters. We derived the masses of the clusters using 
the age-dependent ma s s-to-light ratios of the GALEV mod- 
els dSchulz et all 120021: lAnders & Fritze-v. Alvenslebenl Eo03) 



with Z = 0.004 for the SMC and Z = 0.008 for the LMC. 
For NGC 5236 (M83) we used the ages and masses de- 
rived by Mora et al. (2008, in prep), who kindly provided 
us with a table containing ages and masses of 219 clusters. 
For M33 we used a recent catalogue of 201 clusters pub- 
lished by ISaraiedini & Mancond d2007t) . The M51 data were 
taken from lBastian et al.l (120051) . We deri ve the ages and masse s 
of Antennae clusters from Fig. 2 of Zhang & Falll (1 1999b . 
and converted luminosities to masses using the Bruzual & 
Chariot (1996 , unpu blished) SSP models, which were used by 
Zhang & Fall (1999) to derive the ages. 



consistent with the size of sample prediction without disrup- 
tion (Eqs. I7MT21). The slopes are not exactly +1, but we can 
constrain /mid using Eq.[TT]to /mid < 0.2 for the solar neigh- 
bourhood^ SMC and LMC and f M1D < 0.4 for M33 and M83. 
This rules out the long term 90% ( /mid =0.9) MID proposed by 
Whitmore et al. (2007) as a universal cluster disruption mech- 
anism for these galaxies. For this to be true, M max had to be 
constant with log(age) (Eq. \13\l for all galaxies. 

The trends for M5 1 and the Antennae are clearly different 
from the other five galaxies. The slow or lack of increase of 
log(M max ) for M5 1 in the first ~ 10-20 Myr is consistent with 
the 70% infant mortalit y, independent of m ass, over that time 
range as determined by Bastian et al. (2005). After log(age)=7 
the value of M max is increasing again, confirming that the in- 
fant mortality phase lasts for about 10 Myr only. The fiat slope 
beyond 100 Myr in M51 was inte rpreted as a truncatio n of the 
CIMF around M up =* 5 x 10 s M (bieles et al.ll2006allbl) . Note 
that these three phases can not be derived from the trend of 
M max with log(age) only, since we have only five data points. 
The trend does support the results derive d from the dN / dt dis- 
tribu tion and the luminosity function by Bastian et alJ ((2005) 
and lGieles et alJ d2006albl) 

The fiat relation for the Antennae galaxies is consistent with 
90% disruption each age dex (/mid =0.9) during a Gyr. Note, 
however, that if this r esult is attributed entirely to /mid=0.9 
(Whi tmore et alJ 120071) . which implies that there is no trun- 
cation of the CIMF and that the CFR has been constant ove r 
this age range, which was proposed by I Whitmore et al]d2007l) . 
then based on size-of-sample effects we would expect that the 
galaxy was producing clusters with masses up to 10 9 M Q in the 
oldest log(age) bin, but they have been destroyed due to MID. 
This seems unlikely given what we know about other galac- 
tic mergers. Namely, the major burst of star formation in the 
Antennae s till has to happe n,which will be when the nuclei co- 
alesce (e.g. Cox et alJ2006l) and there are no star cluster known 
with masses in excess of 10 8 M B . 



3.2. Observed trends of log(M max ) vs. log (age) 

In Fig. |2] we show M max as a function of log(age) in bins of 0.5 
dex as open circles. We also show the third most massive clus- 
ter in each age bin (M max ,3 r d) as filled circles, which is expected 
to follow the same relations as the most massive, though with 
less scatter. We note that the trends for the second, fourth and 
fifth most massive cluster are all similar to the trend of M maXi 3 r d. 
The full lines represent the result of fits of log(M max 3,d) with 
log(age) over the first 100 Myr and the dashed lines consider 
a time range of 3 Gyr. The scatter of the data points around 
the fit is most likely due to stochastic effects and should not 
be interpreted as variations in the CFR. The probability density 
function (PDF) of the most massive object drawn from a power- 
law distribution has an asymmetric spread around the mean, but 
in logarithmic units the shape and w i dth PDF is independent 
of the value of M max dLarsenl l2002t iMaschberger & Kroupa 
120071) . Stochastic effects will therefore not affect our results. 

For the solar neighbourhood, SMC, LMC, M33 and M83 
the observed increase of log(M maXj 3 r d) in the first 100 Myr is 



3.3. Deriving the formation rate in bound clusters 

Under the assumption of a constant formation rate, a power-law 
CIMF with index -2 and no disruption of clusters, the values 
of M max are indicative of the amount of stellar mass formed in 
bound clusters. Here we will use the SMC as an example, since 
it has been shown that fo r the clusters in this galaxy the dis- 
ruption time-scale is jong dHodgel 19871 Boutloukos & Lamers 



120031: iGieles et al.l 120071) . In addition, this is one of the only 



galaxies for which a global star formation rate (SFR) has been 
determined and the full galaxy has been imaged for its cluster 
population. This is necessary since M max scales with the total 
number of observed clusters and therefore with the fraction of 
the galaxy that has been imaged. 

For a — 2 the total mass formed in clusters, M tot , can be 
found from 



rM up 
JM min 



M0(M) dM 



(14) 
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(D<1 kpc) 




°; slope= 0.92 
; slope= 0.58 




3 '°; slope= 0.93 
t<10 9 ' 5 ; slope= 0.75 



LMC 




t< 1 a -°; slope= 0.96 
t< 1 9 - 5 ; slope= 0.63 



M33 




t< 1 B -°; slope= 0.76 
t<10 9 ' 5 ; slope= 0.68 
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01 ~z 
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M83 




t< 1 B -°; slope= 0.81 
t< 1 9 - 5 ; slope= 0.83 




t<10 8 -°; slope= 0.45 J 
t<10 9 - 5 ; slope= 0.23 



Antennae 



t<10 8 '°; slope= 0.01 
t<10 9 ' 5 ; slope=-0.21 



6 7 8 9 
log(t/yr) 



10 11 6 



7 8 9 10 
log(t/yr) 



7 8 9 10 
log(t/yr) 



Fig. 2. Evolution of M max with equal size log(age) bins for the seven galaxies in our sample (open circles). The third most massive 
clusters (M maXi 3 r d) in each bin are shown as filled circles. Fits to log(M maXj 3 r d) vs. log(age) on two different age ranges are shown 
as dashed and full lines. 



M m 



\QM n 



(15) 
(16) 



where we again assume M up » M max and M up /M m j n = 10 5 , 
e.g. M min = 10 2 M o andM up = 10 7 M o .Fora = 2 the predicted 



M, 



max, 3rd 



is simply one third of M n 



When representing the data as in Fig. [2] we need to assume 
a CFR, i.e. the amount of mass formed in bound clusters per 
unit of time. With this and the width of the age bins we can 
calculate M tot =CFRx At . We find good agreement between the 
SMC data in Fig.[2]and a prediction for M max and M max ,3rd with 
CFR^ 2 x 10" 3 M o yr _1 . Due to the scatter in the data this is 
only a rough prediction, accurate to within a factor of two. 

We can compare the prediction of t he CFR to the star for- 
mation rate (SFR). lHarris & Zaritskvl (120041) derive a mean 
global star formati on rate of 0. 1 M^ yr -1 based on the field star 
population and IWilke et all d2004h derive SFR=0.05 A^yr" 1 
based on far-infrared observations of the SMC. Our derived 
CFR of 2 x 10~ 3 M yr 1 represents only 2 - 4 % of the total 
SFR. This implies that ~97% of the star formation occurs in a 
dispersed fashion, or, that the infant mortality rate in the SMC 
is close to 97%. This data does not allow to tell these two sce- 
narios apart. 

For the solar neighbour hood, a similar rat io of CFR/SFR^ 
0.05 was found ( Earners & Gielesl2006ll2007l) from a compari- 
son of the star formation rate in (optical detected) clusters rela- 
tive to the star formation in field stars, consistent with the infant 
mortality rate derived by LL03. 



4. Discussion 

4.1. Is there no infant mortality of clusters in some 
galaxies? 

We conclude that the slopes for the first lOOMyr in the solar 
neighbourhood, SMC, LMC, M33 and M83 are consistent with 
a complete lack of mass-independent disruption. This implies 
that the dN/dt of mass limited cluster samples should be fiat 
over this age range and that a luminosity limited cluster sample 
should follow the decline predicted by the fading of clusters, 
i.e. dN/dt oc t~ ( , with - 1.0 < £ < -0.7 dBoutloukos & Lamers 
20031: iGieles et al.l2007l). For the SM C thi s was recently shown 



to be t he case by IGieles et al. I J2007I) and lde Grijs & Goodwin! 
(|2008h . The dN/dt of the LMC clu sters also follows the fad - 
ing prediction in the f i rst lO OMyr dde Griis & Anders! [2Q06h . 



Parmentier & de Griisl (|2008) show that the dN/dt of massive 
clusters is nearly fiat. Both these findings and the near linear 
increase of log(M max ) with log(age) we see in Fig. |2] suggest 
that there can not be 90% MID of clusters in the LMC. 

The dN/dt of M33 c lusters declines as r u±(u 
dSarajedini & M ancone 2007 ). which was interpreted by 
the authors as rapid dissolution. However, the authors also note 
that their sample is luminosity limited. When the sample is 
limited by a detection in a blue filter, such as B or U, the dN/dt 
distribution of a sample that is not affected by any disruption 
declines a s t~ l . This means that both the dN/dt distribution as 
shown by Sarajedini & Mancone (2007) and the nearly linear 
increase of log(M max ) with log(age) (Fig. |2]l are in agreement 
and imply no mass independent disruption (i.e. /mid=0) in the 
first lOOMyr in M33. 
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The flat AN I At distributions of mass-limited cluster sam- 
ples and the results presented in this work are not to be inter- 
preted as no infant mortality of clusters in these galaxies. This 
is because all of the cluster populations used in this study were 
optically selected, i.e. no embedded clusters were included. In 
fact, our findings are in perfect agreement with the 90% infant 
mortality scenario for the solar neighbourhood by LL03, since 
they derived this from a comparison between embedded clus- 
ters to open clusters (see § 13.31 1. Our results suggest that for 
some galaxies the infant mortality of clusters can not be ob- 
served from optically selected clusters. The cluster samples of 
M5 1 and the Antennae are the only ones for which it has been 
shown that there is a steep drop in AN I At around -10 Myr in a 
mass limited, optically selected sample. 

4.2. What happens after 100 Myr? 

The slopes determined over 3 Gyr (Fig. [2]i are shallower for 
most samples, but it is difficult, if not impossible, to come up 
with an explanation where MID disruption starts after 100 Myr. 
There are several culprits for the observed flattening: stel- 
lar evolution combined with "standard" mass-dependent dis- 
ruption, a CIMF that is steeper or truncated at high masses 
(~ 10 6 M ) or a non-constant CFR. Below we will discuss each 
of these effects in turn. 



4.2.1 . The effect of a non-constant cluster formation 
history 

So far we have assumed, for simplicity, that the star/cluster for- 
mation rate has been constant in time. For the age range of 
10-100 Myr this is probably a reasonable assumption for most 
of the galaxies considered in this study. However, for the age 
range of 3 Gyr, as considered in Fig. [2] this assumption may 
not hold and can lead to a confusion between cluster disruption 
and variations in the formation history. For example, for the 
SMC it is generally believed that the assumption of a constant 
CFR is not too bad. Howe ver, when approximating the result 
of lHarris & Zaritskvl d2004 for the global SFR(f) as f, we find 
j] = -0.15. Assuming that the CFR follows the SFR, this par- 
tially explains the small difference between the observed slope 
of +0.75 (Pig. [2j and the predicted slope of +1 (Eq. O, since 
including the effect of this increasing CFR(f) predicts a slope 
of +0.85 (Eq. [8}. 

For the Antennae galaxies the assumption of a constan t 



CFR over the past Gvr (Fall et al]|2005t IWhitmore et alj|200 



is less likely. According to the numerical models of Barne 
dl988l) the first encounter between NGC 4038 and NGC4039 
was around ~ 200 Myr ago and the models of iMihos et al 



( 1993) predicted that the SFR has been increasing by at least 
a factor of five since then. The r l disruption model introduced 



by iFall et al.l d2005l) was derived from only 3 or 4 histogram 
points of the AN I At of the Antennae clusters. A combination 
of 90% infant mortality in the first 10-30 Myr and a CFR that 
started increasing a few 100 Myrs ago could also very well ex- 
plai n the trend in Fi g.[2]fl«fl? the steep decline in AN I At reported 
bv lFall etalJd2005h . 



There will always be a natural bias to studies of star clus- 
ter populations in galaxies which are actively forming stars and 
clusters at the present day, such as the Antennae galaxies, sim- 
ply because they are more attractive to study. The degeneracy in 
this method between changes in the CFR and disruption and/or 
a truncation in the CIMF is hard to quantify. A note of caution 
should be placed that the assumption of a constant CFR can 
lead to false detections of MID. This is not only the case when 
studying AN /At, but as well as M max it vs. log(age). 

4.2.2. Mass dependent disruption 

To flatten the relation between log(M max ) and log(age) by mass 
dependent disruption, the dependence of the disruption time 
(?di s ) on mass has to be weaker than linear, i.e. t&s oc M y where 
y < 1. When fdis K M and M max scales also linear with age 
due to size of sample effects, than all clusters lose the same 
fraction of their initial mass, which would not flatten the rela- 
tion. From obs ervations and theory it is expected t hat f^s scales 
as M 6 dLamers et alj |2005b; Gie les et al. 



2006 cj), from which 



Giele s et alj ([2006a) 



a flattening of M max with age is expected, 
showed that the trends of M max with age for the SMC and LMC 
can be explained by mass loss at old ages due to stellar evolu- 
tion and tidal effects. 



4.2.3. The cluster initial mass function 

Since the measured slope in Fig.|2]is a combination of the index 
of the CIMF and disruption (Eq.[TTIi. the shallower slopes could 
also be explained by a steeper CIMF. 

We note that a slope of +0-75 in Fig. [2] follows from 
a - 2.35 (Eq. 171 iHunter et aT 2003 ). The same high value, 
i.e. larger than 2, for the index of the CIMF was found from the 
analyses of young (< 10 Myr) clusters through the increase of 
M mm as a function o f number of clusters in different galaxies 
( Wei dner et al.ll2004l) . A value of -2.3 was also found for the 
index of the cluster luminosity function through the increase 
of the most luminous cluster as a function of the number of 
clusters dLarsenll2002llWhitmorell2003l:lGieles et alj|2006al) . 

Since these studies consider young clusters, i.e. not affected 
by disruption yet, a Schechter type upper mass limit to the 
CIMF is a logical scenario that can also explain the slopes 
slightly shallower than +1 in Fig. [2] for the fits on the 3 Gyr 
age range. Random sampling from a CIMF that is a power-law 



with a n exponential cut-off around M» = 10 6 M ( Giel es et al 
2006allbh can result in a slope of ~ 0.75 in the log(M max ) vs 



age relation. When sampling cluster masses from such a CIMF, 
the relation of M max with N still resembles a straight line, i.e. 
the cut-off is not detectable as such as long as M max < M„. In 
this scenario slightly steeper slopes at young ages are expected, 
since there the masses are well below the cut-off mass. This can 
also explain the nearly flat slopes of M5 1 and the Antennae, 
since the mass of the most massive cluster is already close to 
M» at youn g ages and can therefore not increase more with age. 
We refer to Gieles (2008) for a more detailed discussion on how 
a physical maximum to the cluster mass can be derived using 
the method presented in § [2] 
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5. Conclusions 

We have studied the evolution of the maximum cluster mass 
{M mdx ) in bins constant in log(age) (i.e. increasing width in lin- 
ear age) in different galaxies. Under the assumption of a con- 
stant cluster formation rate, a power-law cluster initial mass 
function (CIMF) with index -2, and no mass loss of clus- 
ters due to disruptive effects, we predict a linear increase of 
log(M max ) with log(age), a slope of +1, when sampling the 
cluster masses stochastically from the CIMF. Including the ef- 
fects of mass independent disruption (MID) causes the slope to 
decrease in proportion to the fraction of clusters removed per 
age dex, /mid (Eq. 12). This results in a flat relation between 
lo g(M max ) w i th log ( age) when fMip = 0- 9, th e value reported 
by Fa 



(2007). 



1 et al. (2005), Chandar et al. (2006) and Whitmore et al. 



Based on a comparison with observed cluster populations 
in seven galaxies we conclude that: 

1. The predicted linear increase of log(M max ) with log(age) 
for a power-law CIMF with index -2 and no disruption 
is observed in the first lOOMyr for clusters in the so- 
lar neighbourhood, SMC, LMC, M33 and M83. This im- 
plies that no significant disruption has taken place (after 
the clusters have left the embedded phase) and with this 
we can rule out the scenario in which 90% of the clus- 
ters are destroyed each age dex, independent of the cluster 
mass. Comparing the observed slopes with those predicted 
for different mass independent disruption fractions, we can 
constrain /mid < 0.2 for the solar neighbourhood, SMC and 
LMC and / MID < 0.4 for M33 and M83. 

2. For M51 the observed slope i s slightly sha l lower con- 
sistent with the observations of iBastian et alj (|2005) who 
found that 70% of the clusters in M51 do not survive the 
first 10-20 Myr. If we only look at ages between 20 
and 100 Myr the slope is steeper which agrees well with 
no mass-independent disruption, providing additional ev- 
idence that this mass independent disruption phase lasts 
only 10-20 Myr and is due to the expulsion of gas. 

3. In the Antennae galaxies log(M max ) is flat in the first 
100 Myr and even decreasing for the full age range. These 
observations are consistent with the scenario of 90% MID 
each age dex. We note, however, that a non-constant CFR 
or a truncation of the cluster mass function at the high mass 
end can produce the same effect. 

4. The SMC clusters show the best agreement with a pre- 
diction without disruption when assuming a CFR of 2 X 
10 _3 M o yr _1 . This is only 2 or 4% of the global star for- 
mation rate in this galaxy, suggesting that the SMC is ex- 
tremely inefficient in producing bound star clusters, but 
once they are formed they are very stable against disrup- 
tion. 
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